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The construction of a reliable potential for Ge02 , from first-principles, is described. The obtained 
potential, which includes dipole polarization effects, is able to reproduce all the studied properties 
(structural, dynamical and vibrational) to a high degree of precision with a single set of parameters. 
In particular, the infrared spectrum was obtained with the expression proposed for the dielectric 
function of polarizable ionic solutions by Weis et al. [J.M. Caillol, D. Levesque and J.J. Weis, J. 
Chem. Phys., 91, 5544 (1989)]. The agreement with the experimental spectrum is very good, with 
three main bands that are associated to tetrahedral modes of the Ge02 network. Finally, we give a 
comparison with a simpler pair-additive potential. 
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I. INTRODUCTION 

In the vitreous and liquid states, at ambient pres- 
sure, germania (Ge02) (a close structural analog of 
silica (Si02)) forms a tetrahedrally coordinated three- 
dimensional network- . Because of its lower abundance, 
its usage in practical applications is much less widespread 
than for the Si02. Still, Ge02 is used in several fields, 
mainly related to optical technologies. For example, a 
mixture of Si02 and Ge02 allows precise control of re- 
fractive index in optical fibres and waveguides. It is of 
interest to develop simulation methods to allow a detailed 
examination of the local structure in such mixtures and 
to predict the infrared spectrum, since this determines 
the long wavelength limit for their use as optical fibres. 

The essential similarity between the structures of 
glassy Ge02 , Si02 and BeF2 was made clear by the study 
of vibrational properties by Galeener and co-workers^ 
who showed that the inelastic neutron, infrared and Ra- 
man spectra of the different materials were closely re- 
lated. More recent infrared studies of Ge02 have also 
been reported2i>^. Recently a full set of partial struc- 
ture factors were determined in this system by using the 
method of isotopic substitution in neutron diffraction 
experiments^i^ . The results show that the tetrahedral 
network structure is based on corner sharing Ge(Oi/2)4 
tetrahedra with a Ge-0 average distance of 1.73 A and 
with a mean inter-tetrahedral Ge-O-Ge angle of 132°. 
They also show that the topological and chemical order- 
ing in the network display two characteristic length scales 
at distances greater than the nearest neighbour. 

Several molecular dynamics (MD) simulations have 
been undertaken to study the structural and vibrational 
properties of the disordered phases of Ge02. Classi- 
cal MD simulations were performed on glassy and liquid 
Ge02^iS,, using pairwise additive potentials with partial 
charges developed by Oeffner and Elliott (OE) for mod- 
eling the a-quartz and rutilelike phases of Ge02^. More 
recently, first-principles molecular dynamics (FPMD) of 
glassy Ge02 have also been reporte d^°di . In principle, 
the amount of empirical information needed to set up a 



first-principles calculation is minimal and it would nor- 
mally be the method of choice to study the physico- 
chemical properties of condensed phase systems. How- 
ever, it is computationally very expensive compared to 
classical molecular dynamics, which is a major drawback 
when dealing with glassy systems. Classical and FPMD 
simulations were also compared one with each other in 
a combined studyl^ for temperatures above 2530 K. In 
this study apparent inconsistencies between the proper- 
ties predicted by the classical simulations of different au- 
thors are noted. In reality, two different parameter sets 
were proposed by OE: an original one, which was fitted 
from an ab initio energy surface, and a so-called rescaled 
one, which was developed from the previous one in order 
to give a better reproduction of the vibrational proper- 
ties. For example, the partial charge of the germanium 
ion was shifted from 1.5 e in the original potential to 
0.94174 e in the rescaled potential. The apparent incon- 
sistencies arise because different classical potentials were 
being used. Our objective is a single potential which will 
give a good description of structure and dynamics and 
allow long simulations on systems of many atoms. 

The structure of the simulated system can directly be 
compared with the neutron / X-ray diffraction exper- 
iments. Both classical and FPMD simulation provide 
structure factors which are in rather good agreement with 
experiments. For the classical molecular dynamics, the 
most noticeable difference is a shift to high g-values of 
the first sharp diffraction peak (FSDP), which was ob- 
served with both versions of the OE potential^'^. As 
this peak is associated to medium range distances^, this 
means that the Ge-O-Ge bond angle is too large and 
the topology of the tetrahedral network is not correctly 
described. The FPMD simulations provide a diffraction 
pattern with improved agreemenl^, but the small size 
/ timescales explored lead to large statistical uncertain- 
ties, mainly in the region of this FSDP. The vibrational 
properties can also be straightforwardly calculated. The 
vibrational density of states (VDOS) was estimated from 
FPMD and classical MD using the rescaled OE potential, 
with an overall reasonable agreement with the available 
experimental data extracted from inelastic neutron scat- 
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tering experiments. The infrared absorption spectrum, 
which originates from the polarization fluctuations as- 
sociated with motion of the atomic charges, has only 
been reported from the FPMD simulation. The agree- 
ment with experiment^ was overall good, despite a small 
shift of the peaks towards lower frequencies. 

The dynamical properties which involve structural re- 
arrangements are out of reach of the FPMD. The only 
studies on Ge02 employ classical MD, and were aimed 
at determining the diffusion coefficients in the liquid 
stateiii^. Hawlitzky et alJ^ observed important differ- 
ences between their values and those obtained by Mi- 
coulaut et al^; their hypothesis for rationalizing this dis- 
crepancy was inadequate equilibration in the latter. In 
fact, as noted above, it is due to the use of the two differ- 
ent OE potentials in these studies. Indeed, in simulations 
performed with the rescaled OE potential the fluidity of 
the system is enhanced by more than one order of magni- 
tude compared to the results obtained with the original 
OE potential. 

In this work we describe the parameterization of a 
new classical MD interaction potential for Ge02, by us- 
ing some techniques we have developed for halides and 
oxidesi^ii^, which is entirely based on first-principles elec- 
tronic structure calculations. In order to provide an ac- 
curate, transferable description of the interactions, these 
potentials include dipole polarization effects and the ions 
carry full valence charges. The various structural, dy- 
namical and vibrational properties are then compared 
with the corresponding experimental results. In future 
work we will use such potentials to study mixtures and 
the effect of pressure on glassy Ge02 



II. INTERACTION POTENTIAL 
DEVELOPMENT 



The interaction model used in this work is related to 
that used in previous workJ-i^^, with a couple of minor 
differences. The model (known as DIPPIM) includes a 
pair potential, together with an account of the polariza- 
tion effects which result from the induction of dipoles 
on the 0^~ ions. The parameters for these potentials 
were obtained by force-matching them to first-principles 
reference data^. Such an approach was successfully ap- 
plied in the case of other oxide materials^. Tangney 
and Scandoloii also used a similar interaction model to 
study silica, in which all the parameters including par- 
tial charges on the ions were fitted against first-principles 
data. Here we use formal ionic charges (Ge'*"'", O^^) 
which should ensure a better transferability. In the next 
paragraphs we will give a brief description of the model 
used, the first-principles reference calculations and the 
force-matching procedure. 



A. The DIPole-Polarizable Ion Model (DIPPIM) 

The interatomic potential is constructed from four 
components: charge-charge, dispersion, overlap repulsion 
and polarization. The first three components are purely 
pairwise additive: 
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where qi is the formal charge on ion i. The dispersion 
interactions include dipole-dipole and dipole-quadrupole 
terms 
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Here Cg^ and Cg^ 



are the dipole-dipole and dipole- 
quadrupole dispersion coefficients, respectively, and the 
are the Tang-Tonnies dispersion damping function, 
which describe short-range corrections to the asymptotic 
dispersion term. The short range repulsive terms are 
approximately exponential in the region of physical inte- 
rionic separations. The full expression used here for the 
short range repulsion is: 



yrop 



E 



(3) 



where the second term is a Gaussian which acts as a 
steep repulsive wall and accounts for the anion hard core; 
these extra terms are used in cases where the ions are 
strongly polarized to avoid instability problems at very 
small anion-cation separations^^. The polarization part 
of the potential incorporates dipolar effects only. This 
reads: 
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Here ai is the polarizability of ion i, are the dipoles 
and T^^', T^^^ are the charge-dipole and dipole-dipole 
interaction tensors: 



,(5) 

The instantaneous values of these moments are obtained 
by minimization of this expression with respect to the 
dipoles of all ions at each MD timestep. This ensures 
that we regain the condition that the dipole induced by 
an electrical field E is aE and that the dipole values are 
mutually consistent. The short-range induction effects on 
the dipoles are taken into account by the Tang-Toennies 
damping functions: 
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The parameters (f^ determine the range at which the 
overlap of the charge densities affects the induced dipoles, 
the parameters determine the strength of the ion re- 
sponse to this effect. It is important to notice that anion- 
anion damping terms have been taken into account, con- 
trary to what was done in^^. The addition of anion-anion 
damping terms was found to greatly improve the ability 
to match the first-principles data. 



B. The DFT reference calculations 

The parameters in the interaction potential are deter- 
mined by matching the dipoles and forces on the ions 
calculated from first-principles on condensed phase ionic 
configuration o^^'^'^ . Starting from the empirical pair po- 
tential in-, we obtained atomic configurations for Ge02 
by running short MD simulations on small cells (150 
ions); a total of two liquid configurations were obtained 
for Ge02; for each of these, the Hellman-Feynman forces 
acting on individual ions of the simulation cell were cal- 
culated using planewave-DFT code CPMD^^. In all the 
calculations we used norm-conserving pseudopotentials 
and planewave energy cut-offs of 1360 eV; all calcula- 
tions were performed using the generalized gradient ap- 
proximation (GGA) according to the Perdew, Burke and 
Ernzernhof (PBE) scheme. For the calculation of first- 
principles dipoles, the Kohn-Sham orbitals are localized 
via a Wannier transformation to construct maximally lo- 
calized Wannier functions (MLWF) and the dipoles de- 
termined from the positions of the centres of the Wannier 
functions associated with each ior^. 



C. Fitting procedure 



gen terms and rescaled these by the Ge polarizability for 
the anion-cation terms. The Gaussian parameters too 
were added after the fit and then the comparison with 
the first-principles data was run again to check that its 
quality remained unchanged. 

In Figure [T] we report the agreement between the DFT 
calculated forces and dipoles and those predicted by the 
fitted potential, for a set of representative ions; the ab- 
scissa shows the x component of the forces and dipoles, 
while the ordinate is simply the ion number. The quality 
of the representation is quite good and comparable with 
the one obtained in22,; this can be regarded as a very good 
result, considering the simplicity of the model (both ion 
shape deformation effects and quadrupoles are not taken 
into account in the present mode l^^'^^ ). The parameters 
obtained for the DIPPIM potential are summarized in 
Table H 
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The potential parameters are optimized by fitting the 
forces and dipoles, obtained with the DIPPIM poten- 
tials for the reference configurations, to the respective 
results from the DFT calculations; the 2 configurations 
provide a total of about 1800 data points to fit, compris- 
ing three Cartesian force components of each individual 
ion and three components for the dipoles. While most 
of the potential parameters are left as free parameters 
in the fits, there are some exceptions. The polar- 
izability, for instance, was fixed to ao^- =11 a.u., i.e. 
the value Salmon et al obtained experimentally in^. This 
value for the oxide polarizability is also compatible with 
the range of values obtained for this quantity in several 
magmatic melts from first-principles calculations^*'. One 
problem in DFT calculations is the uncontrolled repre- 
sentation of the dispersion interaction. Although disper- 
sion only contributes a tiny fraction to the total energy, 
it has a considerable influence on phase transition pres- 
sures and on the material density and stress tensor. For 
this reason, we decided not to include the dispersion pa- 
rameters in the fits as discussed inH and to add them 
afterwards. We used the parameters in^ for the oxy- 



FIG. 1: Agreement between the DFT calculated dipoles (top 
panel) and forces (bottom panel) and those predicted by a 
fitted DIPPIM potential model for a set of 50 representative 
ions. In each panel the ordinate gives the ion number and the 
abscissa gives the DFT values of the x component of the forces 
and dipoles as points and the fitted values as the continuous 
line. 



III. SIMULATION DETAILS 

We performed MD simulations in the NVT ensem- 
ble through the use of the Nose-Hoover chain thermostat 
method^^i^i. The simulation cell contained 288 oxide and 
144 germanium ions, and its volume was set to 6856 
in order to match the experimental density at 300 K, 
which is 3.66 g.cm^^. We used a time step of 1 fs to in- 
tegrate the equations of motion and the minimization of 
the polarization energy was carried out with a conjugate 
gradient method. The system was studied in both the 
liquid and amorphous phases. 
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TABLE I: Parameters in the repulsive, polarization and dis- 
persive parts of the DIPPIM potential. All values are in 
atomic units. The oxide ion polarizability was set to 11 au. 

0-0 O-Ge Ge-Ge 



A'^ 17.78 40.89 

a'^ 3.893 1.186 

B'^ 50000 50000 

b^^ 1.1 1.75 



CI' 



44 
853 
1.0 
1.0 
1.0 
1.0 



4 
50 
1.0 
1.0 
1.5 
1.5 



2958 
5.96 


/ 




/ 
/ 
/ 
/ 



2-208 



1.977 / 
" 1.770 1.709 / 



For the amorphous phase, the system was first equi- 
librated for 2 ns at the temperature of 4500 K. We 
have then cooled down the system with a cooling rate of 
3.33x10^^ K.s^^, by rescaling the velocities and decreas- 
ing the target temperature of the thermostat in order to 
reduce the temperature by 50 K every 150000 MD steps. 
Finally, a 5 ns long simulation was undertaken at 300 K, 
from which we computed all the data presented in this 
article. This procedure was performed with two differ- 
ent potentials, the DIPPIM, and the original OE, which 
has extensively been used in the literature^. A simula- 
tion at 300 K was also performed with the rescaled OE 
potential (with one of the configurations obtained from 
simulations with the original OE potential as a starting 
configuration) for the study of the vibrational properties. 

For the liquid, a series of long simulations with 
the DIPPIM potential at the temperatures of 3600 K, 
3800 K, 4000 K, 4200 K, 4400 K, 4500 K, 4600 K, 4800 K 
and 5000 K was undertaken. In all the cases, an equili- 
bration run was performed so that the slower species, 
Ge, moved on average, at least a distance of 5.5A; a 
subsequent run of the same length was made to accu- 
mulate enough statistics for the mean-squared displace- 
ments curves. At these temperatures, the runs were be- 
tween 100,000 and 1,000,000 steps. Both the original and 
rescaled OE potentials were also used within the same 
simulation conditions at a temperature of 3600 K. 



IV. STATIC STRUCTURE 

Neutron diffraction experiments with the isotope sub- 
stitution technique allowed Salmon to determine very 
precisely the full set of partial radial distribution func- 
tions (RDF)^. These provide a good check of the validity 



of the potential. 

The experimental RDFs are given in figure [U together 
with the ones obtained from our simulations at 300 K 
involving the DIPPIM and the original OE potential; 
it can be readily observed that our DIPPIM potential 
gives a closer match with the experimental data than 
the OE potential, for which all the characteristic first- 
neighbour distances are slightly overestimated. Concern- 
ing the shape of the function, the peaks obtained for the 
Ge-Ge and Ge-0 partials are too sharp but it has be 
be remembered that other factors, such as system size 
and cooling rate, also play a role in this aspect of the 
comparison. The 0-0 partial RDF is almost in perfect 
agreement with the experimental one, this is very impor- 
tant since the 0^~ anions constitute 2/3 of the atoms of 
the system and their packing arrangement determines the 
arrangement of the tetrahedral network. The OE poten- 
tial gives a reasonable description of the static structure 
(though significantly poorer than the DIPPIM one); one 
might therefore be tempted to use this potential, due to 
its simpler form and faster computational time; we will 
see in the next sections that a reasonable reproduction of 
the structure might not be sufficient for it to be capable 
of predicting the dynamical and vibrational properties. 



V. DYNAMICS 

In this section we will calculate diffusion coefficients for 
the liquid and compare them both with those obtained 
with the OE potentials and with the experimental values. 
The diffusion coefficient of species a can be obtained from 
the slope of the mean squared displacement at long times, 
i.e. 



t-^oo tot 



where 



{rim - ^^f:{\n{t) 
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(7) 



(8) 



The mean squared displacement curves were evaluated 
from a series of long simulations at temperatures between 
3600 K and 5000 K; these are shown on a log-log scale 
on figure [31 No plateau is observed even at the lowest 
temperature, which shows that the system remains liquid 
in the range of the study. 

The diffusion coefficients were extracted using equa- 
tion [T] these are reported in figure H] and compared with 
those obtained with the OE and rescaled OE potentials. 
Figure m shows many interesting features: first of all, our 
DIPPIM potential gives approximately the same value 
for the diffusion coefficients as the original OE poten- 
tial. Hawlitzky et al^ used this potential to calculate 
the diffusion coefficients and showed that the obtained 
values are in good agreement with the experimental ones 
obtained by converting viscosity data; we can therefore 
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FIG. 2: Radial distribution functions obtained from the 
300 K simulations with the DIPPIM and original OE interac- 
tion potentials. These are compared to the experimental ones 
obtained from neutron diffraction experiments^. 



FIG. 3: Time dependence of the mean squared displacement 
curves for both the oxygen and germanium ions. The tem- 
peratures ranges between 3600 K and 5000 K. 



le-03 



. le-04 



le-05r 



le-0( 



-.D 



Ge 



T ' 1 ' r 

OE rescaled . 




DIPPIM 



2.2 2.4 2.6 

1/T (10"V) 



FIG. 4: Diffusion coefficients obtained with a DIPPIM po- 
tential at different temperatures; data points obtained with 
both the original and rescaled OE potentials at 3600 K are 
also presented. 



conclude that our DIPPIM potential gives a good agree- 
ment with the experimental data as well. A direct com- 
parison with experimental data might be done by using 
the only published datum {Dq — 7 x lO"^'^ cm^s~^) on 
the oxygen diffusion constant^^, at T—IAAO K. However, 
such a low diffusion coefficient would involve an impossi- 
bly long simulation run. 

It is also interesting to see that the rescaled OE poten- 
tial yields too fluid a melt if compared with the original 
OE and DIPPIM potentials and with experimental data. 
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This potential was used by Micoulaut et alF- and they 
indeed obtained diffusion coefficients which were more 
than one order of magnitude larger than the values ob- 
tained by Hawlitzky et ali^^ (who incorrectly explained 
this discrepancy in terms of lack of equilibration of the 
simulations, it is obvious here that the real reason was 
the difference between the two potentials) . The so-called 
rescaled potential was obtained by taking the normal OE 
potential and rescaling the ionic charges, the van der 
Waals coefficient and the repulsive parameters in order 
to get a better agreement with the experimental vibra- 
tional spectrum for the crystal. To do this the parame- 
ters which are responsible for the strength of the the 0-0 
and 0-Ge short range repulsion were lowered by a factor 
of approximately 2.5; as a consequence, these ions can 
diffuse much more easily in a softer environment. 



VI. INFRARED SPECTRUM 

The infrared spectra of ionic melts with polyvalent 
cations exhibit discrete bands attributable to the vibra- 
tional motion of the local coordination complexes around 
the cations. They originate from the polarization fluc- 
tuations associated with motion of the ionic charges. 
The inclusion of polarization effects for the oxide ions in 
our model may influence the predicted spectrum in two 
ways2^. First, the interactions of the oxide ion dipoles 
may alter the local structure of the network and the 
strength of the bonds, which may introduce a shift of 
the vibrational frequencies. Second, the induced dipoles 
will themselves be responsible for absorption, as they too 
contribute to the total polarization fluctuations. 
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FIG. 5: Imaginary part of the dielectric function, calculated 
for simulations carried out with the DIPPIM potential, com- 
pared to FPMD simulationiS' and experimental resulta^i^. 

The absorption coefficient in the presence of these ex- 
tra moments is calculated from the imaginary part of 



TABLE IL Characteristic frequencies of the imaginary part of 
the dielectric function as obtained from the DIPPIM model, 
compared to previous FPMD simulations^"'^^ and to experi- 
mental resulta^^. All the frequencies are given in cm~^. 



DIPPIM 


300 


601 


894 


FPMD^a 


244 498 


787 


Exp. 1^ 


280 


567 858 


Exp. 2^ 


315 


585 


915 


Exp. 3** 




560 


870 



the total dielectric function [n{i>)a{v) = 27rz/9(e(z/))], 
which can be determined following Caillol, Levesque and 
Weis21^ as 



where 



+2(M- J),+ 
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2ttv 



2mv{m ■ M)^ (9) 



(J ■ J). 



J(i) • J(0))di, 



(10) 



J{t) is the charge current J{t) = X]i3=i Qi'^i{t) and M(t) 
is the total system induced dipole moment, M(t) — 

Ef=iM.(t). 

Figure [5] shows the imaginary part of the total di- 
electric function calculated for simulations carried out 
with the DIPPIM potential, compared to the one ob- 
tained from FPMD simulation o^°di and from several 
experiment32i>^ (the spectrum proposed by Kamitsos et 
al. shows the absorption coefficient, here we transformed 
it into the dielectric function to facilitate comparisons). 
Three main bands are observed on the spectra, in good 
agreement with the other studies. The corresponding 
characteristic peak frequencies for the absorption spec- 
tra are summarized in table |TT1 The agreement is very 
satisfactory; the frequencies obtained match very closely 
to the various experimental results (even though these 
show a significant scatter), while the FPMD simulation 
seems to systematically underestimate them. The rela- 
tive intensities of the bands is also in good agreement 
with experiment, though the highest frequency band at 
around 894 cm~^ is not intense enough in our spectrum. 

The importance of polarization effects in determining 
the relative intensities of the bands can be demonstrated 
by separating the various contributions to the absorption 
spectrum. In the case of glassy silica^^ and beryllium 
fluoride'^^, it was observed that the interference between 
the induced dipoles and permanent charge contributions 
to the total polarization, contained in the (M • J),^ cross 
term, is responsible for the changes of the relative intensi- 
ties. Figure E] shows the decomposition of the absorption 
spectrum in the case of glassy Ge02. The two main con- 
tributions are the charge fluctuation and the cross term; 
both show bands at the same intensities. Note that the 
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cross-term strongly reduces the intensity of the two low 
frequency bands relative to that which would be obtained 
from the charge fluctuations alone, whereas for the high- 
frequency band the cancellation is much weaker. For the 
high frequency band the relationship of the charge-charge 
and cross-terms is different to that found previously for 
Si02 and BeF ^^^i^^ . In the latter cases, the charge- 
charge and cross-terms had the same sign, so that the 
net band intensity was slightly larger than that predicted 
by the charge-charge term alone. This difference in the 
behaviour of the calculated spectra might arise from the 
inclusion of an anion-anion damping term in the polar- 
ization part of the interaction potential here, which was 
not the case for the previous studies. 




200 400 600 800 1000 1200 

V (cm ) 



FIG. 6: Decomposition of the total IR absorption spectrum 
in its three component parts. 

We can now attempt to associate these IR frequencies 
with the vibrational modes of the system. The upper part 
of figure[7|shows the vibrational density of states (VDOS) 
of each ion. Some of these features can be associated to 
the vibrations of the Ge04 tetrahedra by comparing them 
with the tetrahedral-VDOS as obtained from the correla- 
tion functions of the velocities of various local tetrahedral 
symmetry coordinates. The latter would correspond to 
the velocities of Ge04 normal modes, if the tetrahedra 
were isolated {i.e. not mechanically linked in a network). 
For a tetrahedral molecule, vibrational normal coordi- 
nates vi {Ai, symmetric stretch), V2 {E, bend), 1^3 {F2, 
asymmetric stretch) and 1^4 {F2, bend) are expected^. 
We may obtain velocities associated with the symmetry 
coordinates of each Ge04 unit. For example, for the sym- 
metric stretching motion 



a=l-4 



(11) 



where i labels the Ge of a tetrahedral complex, and ia 
the four oxide anions at the vertices, vf^ is the projection 



of the relative velocity of ia along the i ^ ia bond, i.e. 

vfa = (v^a - Vi) • {r,a - r^). (12) 

Similar expressions may be written down for the veloc- 
ities of the other symmetry coordinates^. The corre- 
sponding DOS is then obtained by a Fourier transform 
of the corresponding velocity autocorrelation function: 



DOSa(j^) = 5R / 
Jq 



{vfit)-vfmdt. 



(13) 




1200 



V (cm ) 



CO 
O 
Q 
> 



T3 

lU 

H 



— E -bend 

■ ■ ■ Aj -stretch 

- F^-bend 

■ - F^-stretch 




200 400 600 800 1000 1200 
V (cm ) 



FIG. 7: The vibrational densities of states, for the DIP- 
PIM, of both type of ions (top panel) , and their projections in 
terms of the A (symmetric stretch), F2 (asymmetric stretch), 
E (bend), and F2 (bend) modes of the Ge04 tetrahedra (bot- 
tom panel). 

Were the network to vibrate as a collection of isolated 
tetrahedra, the spectra of each tetrahedral-VDOS would 
exhibit a single peak corresponding to the characteris- 
tic frequency of the corresponding normal mode of vi- 
bration. The four tetrahedral-VDOS are shown on the 
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lower part of figure [71 it is clear that it is not the case 
here. Each spectra consists in the superposition of sev- 
eral bands, which shows the existence of some important 
coupling between the symmetry coordinates. It is there- 
fore difficult to assign the IR features to the different 
modes. In particular, all of the three bands observed in 
the IR spectrum appear in the tetrahedral-VDOS of the 
bending E and F2 modes. The other F2 modes, which is 
associated to some kind of stretching of the tetrahedra, 
only exhibits the medium and high frequency modes. Fi- 
nally, the Ai stretching spectrum exhibits two bands at 
frequencies centered on the values of 510 and 970 cm~^, 
which do not appear on the IR spectrum. Such a mode of 
this symmetry is indeed not expected to be IR active for 
an isolated tetrahedral molecule, since it corresponds to 
the symmetric breathing mode. Coming back to the total 
VDOS of each ion, the low frequency vibrations are not 
associated to any of the tetrahedral modes. In fact these 
are due to some slow collective bending of two tetrahedra 
linked by a bridging oxide anioiii^, so that they would not 
appear in the DOS of a single tetrahedron. 

The vibrational properties of the system as obtained 
with the original and rescaled OE interaction potentials 
have also been computed. As these potentials do not 
include any polarization effect, it is only possible to com- 
pare the charge-charge current contributions to the IR 
absorption spectrum with the one obtained with the DIP- 
PIM potential. The obtained spectra are given on fig- 
ure [5] Unsurprisingly, the original OE potential seem to 
fail completely in describing the vibrational properties of 
Ge02; this is the reason why the rescaled OE potential 
had to be developed. For that potential, the agreement 
is better, in particular for both the small and high fre- 
quency bands. Concerning the medium frequency band, 
observed at around 604 cm~^ for the DIPPIM potential, 
it does not appear for the rescaled OE potential. Instead 
there appears a shoulder to the small frequency band, 
for frequencies ranging from 400 to 500 cm~^. As these 
bands were associated to the F2 stretching and bending 
modes on the tetrahedral-VDOS, this means that this 
modes are not correctly depicted by the rescaled OE po- 
tential. 



VII. CONCLUSION 

In conclusion, the obtained potential seem to repro- 
duce all the studied properties (structural, dynamical 
and vibrational) to a high degree of precision. Since 



there was no reference at all to any experimental data in 
the parameterization of this potential, this represents a 
strong test of the model reliability. The study of pressur- 
ized germania is among our top priorities for future work. 
The potential is of a similar form to that introduced to 
describe silicates^; since these potentials have proved to 
be transferable, we will be able, in the future, to focus on 
Ge02-Si02 mixtures as well. These mixtures, as already 
mentioned in the introduction, attract great interest as 
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FIG. 8: Comparison of the charge-charge current contribu- 
tions to the infrared spectrum obtained with the three differ- 
ent potentials (the two OE curves were rescaled with a factor 
go/4, where qo is the partial charge carried by the oxide ion 
in the corresponding potential. 



they are widely used in optical fibers and waveguides. 



Acknowledgements 

We thank Philip Salmon for providing his diffraction 
data. Calculations were carried out at the EaStCHEM 
computing facility (http://www.castchem.ac.uk/rcf). 
DM thanks the EPSRC and the CMPC for his PhD 
funding. This work was carried out under the HPC- 
EUROPA-f-l- project (project number: 211437), with the 
support of the European Community - Research Infras- 
tructure Action of the FP7 "Coordination and support 
action" Programme. 



^ M. Micoulaut, L. Cormier, and G.S. Henderson. The 
structure of amorphous, crystalline and liquid Ge02. J. 
Phys.iCondens. Matter, 18:R753-R784, 2006. 

^ F.L. Galeener, A.J. Leadbetter, and M.W. Stringfellow. 
Comparison of the neutron, raman, and infrared vibra- 
tional spectra of vitreous Si02, Ge02 and BeF2. Phys. 



Rev. B, 27(2) 11052-1078, 1983. 

^ E.l. Kamitsos, Y.D. Yiannopoulos, M.A. Karakassides, 
CD. Chryssikos, and H. Jain. Raman and infrared struc- 
tural investigation of 2;Rb20.(l — 2;)Ge02 glasses. J. Phys. 
Chem., 100(28): 11755-1 1765, 1996. 

^ P.V. Teredesai, D.T. Anderson, N. Hauser, K. Lantzky, 



and J.L. Yarger. Infrared spectroscopy of germanium diox- 
ide (Ge02) glass at high pressure. Phys. Chem. Glasses, 
46(4) :345 349, 2005. 

P.S. Salmon, A.C. Barnes, R.A. Martin, and G.J. Cuello. 
Glass fragility and atomic ordering on the intermediate and 
extended range. Phys. Rev. Lett, 96(23):235502, 2006. 
P.S. Salmon, A.C. Barnes, R.A. Martin, and G.J. Cuello. 

Structure of glassy Ge02. J. Phys.iCondens. Matter, 
19(41):415110, 2007. 

M. Micoulaut, Y. Guissarii, and B. Guillot. Simulated 

structure and thermal properties of glassy and liquid ger- 

mania. Phys. Rev. E, 73(031504), 2006. 

J. Peralta, G. Gutierrez, and J. Rogan. Structural and 

vibrational properties of amorphous Ge02: a molecular 

dynamics study. J. Phys.:Condens. Matter, 20(14):145215, 

2008. 

R.D. Oeffner and S.R. Elliott. Interatomic potential for 
germanium dioxide empirically fitted to an ab initio energy 
surface. Phys. Rev. B, 58(22):14791-14803, 1998. 
L. Giacomazzi, P. Umari, and A. Pasquarollo. Medium- 
range structural properties of vitreous germania obtained 
through first-principles analysis of vibrational spectra. 
Phys. Rev. Lett, 95(7):075505, 2005. 
L. Giacomazzi, P. Umari, and A. PasquareUo. Vibrational 
spectra of vitreous germania from first-principles. Phys. 
Rev. B, 74(15) ;155208, 2006. 

M. Hawlitzky, J. Horbach, S. Ispas, M. Krack, and 
K. Binder. Comparative classical and ab initio molecular 
dynamics study of molten and glassy germanium dioxide. 
J. Phys.-.Condens. Matter, 20(28):285106, 2008. 
P.A. Madden, R. Heaton, A. Aguado, and S. Jahn. 
From first-principles to material properties. J. Mol. 
Struct. -THEOCHEM, 771:9-18, 2006. 

S. Jahn and P.A. Madden. Modeling Earth materials from 
crustal to lower mantle conditions: a transferable set of 
interaction potentials for the CMAS system. Phys. Earth 
Planet Inter., 162:129-139, 2007. 

M.J. Castiglione, M. Wilson, and P.A. Madden. Polar- 
ization effects in the simulation of lead (II) fluoride. J. 
Phys.:Condens. Matter, 11:9009-9024, 1999. 
M.J. Castiglione and P.A. Madden. Fluoride ion disor- 
der and clustering in superionic pbf2. J. Phys.: Condens. 
Matter, 13(44):9963, 2001. 

P. Tangney and S. Scandolo. An ab initio parametrized 
interatomic force field for silica. J. Chem. Phys., 117(19): 
8898-8904, 2002. 

The CPMD consortium. CPMD Version 3.x, 

http://www.cpmd.org. MPI Fiir Festkorperforschung 



and the IBM Zurich Research Laboratory. 
^■^ A. Aguado, L. Bernasconi, S. Jahn, and P.A. Madden. 

Multipoles and interaction potentials in ionic materials 

from planewave-dft calculations. J. Chem. Soc., Faraday 

Discuss., 124:171-184, 2003. 
^° M. Salanne, R. Vuilleumier, P.A. Madden, C. Simon, 

P. Turq, and B. Guillot. Polarizabilities of individual 

molecules and ions in liquids from first principles. J. 

Phys.: Condens. Matter, 20(49) :494207, 2008. 

21 S.T. Norberg, I. Ahmed, S. Hull, D. MarroccheUi, and P.A. 
Madden. J. Phys.: Condens. Matter, to be submitted. 

22 M. Wilson, S. Jahn, and P.A. Madden. The construc- 
tion and application of a fully flexible computer simula- 
tion model for lithium oxide. J. Phys.: Condens. Matter, 
16(27):S2795-S2810, 2004. 

2^ S. Nose. A molecular dynamics method for simulations in 
the canonical ensemble. Mol. Phys., 52:255 268, 1984. 

2" G.J. Martyna, M.L. Klein, and M. Tuckerman. Nose- 
hoover chains: the canonical ensemble via continuous dy- 
namics. J. Chem. Phys., 97(4):2635 2643, 1992. 

•^^ T. Tokuda and W.D. Kingcry. Oxygen diffusion and va- 
porization rates for liquid and solid germanium dioxide. J. 
Appl. Phys., 34:2104-2105, 1963. 

2® M. Wilson, P.A. Madden, M. Hemmati, and C.A. An- 
gell. Polarization effects, network dynamics, and the in- 
frared spectrum of amorphous sio2. Phys. Rev. Lett, 
77(19):4023-4026, 1996. 

J.M. Caillol, D. Levesque, and J.J. Wcis. Electrical prop- 
erties of polarizable ionic solutions. I. Theoretical aspects. 
J. Chem. Phys., 91(9):5544-5554, 1989. 

2* J.M. Caillol, D. Levesque, and J.J. Weis. Electrical proper- 
ties of polarizable ionic solutions. II. Computer simulation 
results. J. Chem. Phys., 91(9):5555-5566, 1989. 

2^ R.J. Heaton, R. Brookes, P.A. Madden, M. Salanne, C. Si- 
mon, and P. Turq. A first-principles description of liquid 
BeF2 and its mixtures with LiF: 1. potential development 
and pure BeFa. J. Phys. Chem. B, 110(23):11454-11460, 
2006. 

^° E.A. Pavlatou, P.A. Madden, and M. Wilson. The inter- 
pretation of vibrational spectra of ionic melts. J. Chem. 
Phys., 107(24): 10446-10457, 1997. 

^1 K. Trachenko and M.T. Dove. Floppy modes in silica glass 
under pressure. J. Phys.:Condens. Matter, 14:1143-1152, 
2002 . 

'^2 S. Jahn and P.A. Madden. Structure and dynamics in liq- 
uid alumina: simulations with an ab initio interaction po- 
tential. J. Non-Cryst Solids, 353(32-40) :3500-3504, 2007. 



